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rH ■ Abstract We propose an analysis technique for the exceptional points (EPs) occur- 

^ . ring in the discrete spectrum of open quantum systems (OQS), using a semi-infinite 

chain coupled to an endpoint impurity as a prototype. We outline our method to lo- 
VO . cate the EPs in OQS, further obtaining an eigenvalue expansion in the vicinity of 

the EPs that gives rise to characteristic exponents. We also report the precise number 
of EPs occurring in an OQS with a continuum described by a quadratic dispersion 
,^ ■ curve. In particular, the number of EPs occurring in a bare discrete Hamiltonian of 

M-i dimension np is given by nuinu — 1); if this discrete Hamiltonian is then coupled 

to continuum (or continua) to form an OQS, the interaction with the continuum gen- 
erally produces an enlarged discrete solution space that includes a greater number of 
r^ ' EPs, specifically 2"'^(ric + nD)[2"'^(nc + tt-d) — 1], in which nc is the number of 

Q^. (non-degenerate) continua to which the discrete sector is attached. Finally, we offer a 

heuristic quantum phase transition analogy for the emergence of the resonance (giv- 
ing rise to irreversibility via exponential decay) in which the decay width plays the 
^ ■ role of the order parameter; the associated critical exponent is then determined by the 

0> , above eigenvalue expansion. 
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1 Introduction 

Exceptional points (EPs) represent defective points occurring in the parameter space 
of some finite-dimensional Hamiltonian Ho at which two (or possibly more) of the 
system eigenvalues coalesce while, simultaneously, the usual self-adjoint property of 
the Hamiltonian fails to hold [Q. EPs serve to characterize the qualitative aspects 
of the spectrum in their vicinity and have been studied in a wide array of physical 
contexts, including experimental studies of the field modes in a microwave cavity JJ] 
13 as well as a chaotic optical microcavity JU . In the context of scattering theory an 
EP at which two eigenvalues coalesce appears as a double pole of the S-matrix; for 
examples see the analysis of nuclear spectra, laser-induced continuum structures in 
atoms and the transmission through quantum dots that is reported in the review Q. 
Further theoretical studies of exceptional points in various physical contexts include 
superconducting wires J61, electron scattering on helium Q, photodissociation of 
H2 El, the hydrogen spectrum under crossed electric and magnetic fields 191, the 
scattering of a beam of particles in a double barrier potential jTOll , and nuclear spectra 
studies relying on random matrix theory ITTl . A quantum phase transition (QPT) 
interpretation has also been offered based on a Coulomb interaction analogy for the 

EPs ESI. 

In physical situations the finite Hamiltonian Ho is usually embedded in a much 
larger system (i.e. open quantum system (OQS)) such that the discrete spectrum of 
Hd is perturbed by the continuous degrees of freedom associated with the larger 
system. The continuum might represent a measuring device, such as the leads used 
to probe the cavity modes in Refs. |l2][3]. Note, however, that the EPs in this case are 
not merely shifted from those of the discrete space by itself; in fact, the interaction 
with the continuum in general gives rise to a greater number of EPs inhabiting an 
effective finite-dimensional Hamiltonian with dimension larger than H^ alone (we 
demonstrate this subtle point for our prototype model below). 

Here we propose a method for the discovery and spectral analysis of EPs in OQS. 
This analysis incorporates a certain eigenvalue expansion that generically occurs in 
the vicinity of any EP 11731 161171. Further, we offer our own QPT analogy that relates to 
the emergence of the resonant state in an open system and hence acts as a mechanism 
for spontaneous time-symmetry breaking in nature, with critical exponents given by 
this eigenvalue expansion. While we will demonstrate this process for a specific ex- 
ample, we emphasize that a similar effect occurs in many more general OQS models. 

The Hamiltonian we have in mind may be generically written as H = Hu + Hq + Hue', 
here Hq gives the continuum (open) component of the system coupled to a discrete 
level or series of levels H^ through the coupling Hamiltonian Hoc- We will focus on 
the case Hu = ej d''d, in which the discrete component is written in terms of a single 
level coupled to the continuum via Hqc = g Vdc with g a dimensionless coupling 
constant. However, we make this specification with the understanding that this single 
level may simply be one particular discrete level, with others absorbed into He, and 
the general ideas from our analysis below still applies in the case of multiple levels 
IT]. We assume that no many -body interactions are present. 

This Hamiltonian (or close generalizations) may be applied in the description 
of a wide range of physical phenomena, for example the interaction of a two-level 
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oscillator with the field modes in a waveguide llT4l or the effect of a quantum dot 
impurity on quantum wires 11511171 and has been studied extensively ll5l [T8l[T9ll20l 
l2ni22l . We comment on the expected general applicability of our method below. 

Defining the Green's function G{z) = [z — H)~ we obtain the discrete spec- 
trum for H as the poles of {d\ G{z) \d), where \d) is the eigenvector of Hu with 
the eigenvalue e^. For the simplest cases, this spectrum appears as the solutions of a 
dispersion equation in the form 

z-ei = S{z). (1) 

with the self-energy function S{z) oc g^. The open system allows for multiple solu- 
tions that may include complex solutions (i.e. resonances) (231. 

We outline our analysis technique in this paper for the case that the dispersion 
equation is given by Eq. ([T). While our method is oriented towards analytically solv- 
able models, it nevertheless should prove useful in some numerical situations assum- 
ing that we may at least explicitly write the self-energy function S{z) in integral 
form. Further, our method should in principle be portable to other OQS such as the 
scattering problem in Ref. ifTOl . For systems involving a more complicated struc- 
ture for Hy) with multiple bound states, one generally obtains a matrix equation ll24ll 
equivalent of Eq. ([TJ, in which case it should again be possible to generalize the 
methods described below. However, for systems involving significantly more com- 
plex interactions between the discrete sector and the fields, such as Ref. IJO], this may 
be difficult to implement in practice, as well for the case involving multi-particle 
interactions. 

Let X generically represent any two Hamiltonian parameters, such as a; = {g, ej} 
above. An EP is a point x in parameter space at which two eigenvalues Zi (x) and 
Zj{x) (though possibly more) coalesce as Zi{x) = Zj{x) = Zc, while their eigenvec- 
tors also collapse into one another. Note that we have defined z^ as the coalescence 
point in the complex energy plane. The coalescence occurs because the two eigenval- 
ues share a common branch cut in the complex parameter space. The branch cut is 
described by the system discriminant D{x), which is obtained in the case of a closed 
system from the determinant of the finite Hamiltonian; in the case of an OQS we may 
obtain D{x) as the determinant of a finite-dimensional effective Hamiltonian i^eff 
that describes the spectrum in the discrete subspace (iJeff may be obtained by apply- 
ing the Feshbach projection operator technique |5 1 or by other methods II23II25I as we 
demonstrate in App. |A]i. Adiabatic encirclement of the defective point x in complex 
parameter space will exchange the two eigenvalues Ifni9l ll3l . 

In a more general context, following Kato's formalism llT3l we may group the 
system eigenvalues in the vicinity of a given EP x into cycles as 

{zi{x),...,Zp^l{x),Zp{x)}, {Zp+l{x),...,Zp+q{x)},... (2) 

Upon a single adiabatic revolution of the system parameters about x the eigenvalues 
will be mapped into one another according to 

{zi{x),...,Zp^i{x),Zp{x)) ->■ {z2{x),...,Zp{x),zi{x)} (3) 

and so forth. Hence we must complete p revolutions in order for this cycle to return 
to its original configuration, where p is the period (the number of elements in a given 
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cycle). Note that we have p = 2, the simplest case, for our prototype model starting 
from Sec. [51 however we maintain that p may take on arbitrary values p >2 for the 
time being. Further, we call the energy value z^ at which the eigenvalues coalesce the 
center of the cycle. 

We offer that the elements of the cycle may be expanded in the vicinity of the EP 
in a fractional power series in the form 

Zh{x) = z, +ai,,, (/j(x))'/P + a2.h ifjix)f + ■■■ (4) 

Here ai^h = Pi exp{2TTih/p) with /3; an unknown coefficient to be determined below 
and h ~ 0,1, . . .p — 1. The functions fj{x) are polynomials obtained by factoring 
the discriminant D{x) such that EPs associated with a shared branch cut have been 
placed in a common factor (an example appears below). Our expression Eq. (|4|i is a 
mild generalization of that appearing in Ref. |[T3l . in which only linear polynomials 
fj {x) = X — X are considered. Note that in the case of an ordinary Taylor expansion 
the distinction would be largely irrelevant, as in that case the expression incorporating 
higher order polynomials could just as well be re-grouped to give the more familiar 
form written in terms of linear factors. However, the fractional powers appearing in 
Eq. (IDi make such a re-grouping impossible in the present case; hence here it is best 
to work with the most general expression possible. 

Having introduced our essential framework for the EPs, we proceed as follows. 
In Sec. 12] we briefly present our technique for locating and analyzing EPs in OQS. 
We turn in Sec. [3] to the specific case of a semi-infinite tight-binding chain with an 
endpoint impurity in order to provide an example where the EP analysis may be 
performed in a simple analytic context with closed-form expressions. In Sec.|4]we 
present a detailed look at the eigenvalue spectrum for the case g < 1/a/2, in which 
a resonance is present in our prototype model for certain values of e^. We further 
obtain the fractional power series Eq. (HI for our prototype for values of ej in the 
vicinity of the point in parameter space at which the resonance appears. We also take 
the opportunity to comment on the number of EPs occurring in OQS, dependent on 
the form of the dispersion in the leads. In Sec.|5]we present our QPT analogy for the 
appearance of the resonance and the consequent breaking of time-reversal symmetry. 
We summarize our work and make final conclusions in Sec. |6l 

We also demonstrate how to obtain an isospectral effective Hamiltonian for our 
OQS in App. lAl and in App.fBjwe comment on the number of solutions appearing 
in an OQS described by a given effective Hamiltonian. The discussion around the 
effective Hamiltonian elucidates the relationship between the typical OQS we have 
in mind and the T'T-symmetric Hamiltonians that have been the subject of numerous 
studies in the last decade or so. 



2 Analysis technique for EPs in OQS 

We now briefly outline our general method for finding the EPs in OQS as well as 
finding the analytic form of the expansion Eq. (JDi. We emphasize that this technique 
should remain useful as a numerical tool in the case of systems significantly more 
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complex than our prototype model introduced in Sec.|4] We comment on some other 
methods for locating the EPs at the end of this section. 

For our approach, we take the derivative of Eq. ([T]) with respect to the impurity 
energy ea and re-arrange to obtain 



dzo dzo/dci' 



(5) 



As the eigenvalue derivative diverges at the EP Q according to dzi/d e^ 

(which can be seen immediately from Eq. (|4]i since p > 1) we may evaluate Eq. (O 

at the EP as 

dS{za) 



dzo 



= 1, (6) 



where Zc is the center associated with the EP l^. This gives us a condition to determine 
the value of the center Zc in terms of the impurity-continuum coupling parameter gii 
We determine the location of the EP in the parameter space (cd, .9) by plugging z^ into 
the dispersion Eq. ([T) and solving for ej for a fixed value of g (or vice versa g for fixed 
Ed); this result naturally occurs as one of the discriminant polynomial factors fj{ed, g) 
in Eq. (HI. While finding the EP, note that a sign ambiguity will arise associated with 
the branch cut; this can be overcome by performing a consistency check on the value 
of Zc obtained in Eq. (|6]l. 

Having obtained z^ along with the relevant EP, we further obtain the higher order 
terms in Eq. (|4]i to arbitrary precision by expanding Eq. (jSj in powers of (/j (ej , 5) ) 
to obtain conditions on the coefficients aiji in terms of the system parameters. The 
remaining unknown is the period p, which is obtained by integrating the logarithm of 
the Green's function G{z) around a closed contour Cg^ in (energy) z-space about Zc 
as 

Note that once the location of the EPs are determined, the value of p could be found 
numerically in a straight-forward manner using this equation. 

Clearly the above method for locating the EPs is only useful in the case that 
we can write some type of closed form expression for S{z). In the case that this is 
not possible, several purely numerical approaches have been proposed in the liter- 
ature II26II27||281 . However, the numerical method presented in Ref. f26\ requires a 
finite matrix as input, while an OQS necessarily contains a continuum element. This 
problem may be overcome by introducing an isospectral finite-dimensional effective 
Hamiltonian H^ff, as in App.|A] This introduces yet another issue; the method in Ref. 
1261 assumes the relevant eigenvalue equation takes the form Htpz = zip^ while the 
effective Hamiltonian method results in a more generalized eigenvalue equation of 
the form Hgf{{z)'i/jz — zip^, as demonstrated in App. lAl Hence one would need to 
check carefully if this method can be appropriately generalized. 



' Note that Eq. ^6) is the equivalent of Eq. (16) in Ref. 1101 : hence the technique described here below 
Eq. (6) could be carried over to the scattering problem as well. 
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3 Prototype model: application to a semi-inflnite chain with end-point impurity 

We now apply the techniques above to a semi-infinite chain with an end-point im- 
purity as an illustration of the emergence of a resonant state, which occurs at a real- 
valued exceptional point. Further on we will offer our QPT analogy for this process, 
with the characteristic exponent for the order parameter (decay width) read off from 
Eq. (|4|i. Our model Hamiltonian (see Fig.fTfa)) takes the form 

1 " 
H = tid^d - 2 X! ('^nC«+i + cl+ic« j - -7= [c\d + d^cij , (8) 



=1 



in which the coefficient of the second term sets the energy scale for this paper. This 
model is exactly solvable and has been the subject of numerous studies Ill5lll6ll20l 
l29l . The site-to-site hopping in the second term of Eq. ^ is easily diagonaUzed by 
introducing a Fourier transform 

which reveals the dispersion Ck = — cos k, covering k e [0, tt] in the continuum 
limit A^ — > oo. The third term gives the impurity coupling to the end-point of the 
chain with strength —g/^/2. We obtain the discrete spectrum for the fully diagonal- 
ized model from the poles of {d\ G{z) \d) = {z — e^ — S{z))~ . Since our model 
is purely quadratic, this expression may be evaluated by a re-summation following a 
straight-forward operator expansion in terms of the coupling. Performing this evalu- 
ation yields the equation for the poles as 



e, = Siz)^g^z-Vz^-l) (10) 



(we present an alternative derivation of this result in App.lAJi. We may square out the 
root in this expression to obtain a quadratic Pq{z), then following Ref. If35]| we obtain 
the system discriminant from the solutions of Pq{z) — P'{z) = as 

D(ed;5) = -4//i(g)/2(6d,.g), (11) 

in which we define the polynomial factors fi{g) = 1 — 2g^ and /2 (ed, .?) = Cd^ ^ (1 ^ 
2g^) mentioned above. The polynomial /2(ed, g) describes a branch cut in complex 
Ed space extending from e_ to e_|- with 



e± 



= ±VT^V, (12) 



revealing Cd = e± as two EPs written in terms of fixed g. We note immediately that 
for g < I/V2 {g > l/\/2) these EPs are real (imaginary); for present purposes we 
will focus on the real case. Meanwhile the polynomial fi{g) reveals that g = l/-\/2 
is (coincidentally) some special point in its own right (this will be further considered 
in a separate publication). 

Before proceeding with a closer inspection of the eigenvalues, we take advantage 
of an opportunity here to comment on the number of EPs appearing in certain kinds 
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Fig. 1 (a) Prototype model geometry: semi-infinite chain coupled with a discrete endpoint impurity; the 
nearest-neighbor coupling is given by —1/2 and the chemical potential of the sites along the chain is 
eo = 0. The impurity has chemical potential ej and is attached to the endpoint site 1 with coupling 
g' = g/\/2- (b) Diagrammatic representation of the finite-dimensional effective Hamiltonian _ffeff (z) 
derived in App.|A] The effect of the chain on the impurity site is captured here by the self-energy function 



of OQS. According to Ref. ifTSl . the number of EPs r appearing in the spectrum of a 
discrete Hamiltonian of dimension n is given by r = n{n—l). However, the situation 
for OQS is a bit more subtle. In the case of our prototype system, we immediately 
notice that the number of EPs is not fixed by the dimensionality of the bare discrete 
sector Hd, according to which there would be zero EPs (after all, a single eigenvalue 
cannot coalesce into itself). To elucidate the situation, it is useful to consider the 
effective Hamiltonian (see Fig.lllb)) 



H,niz) 



9 

\S{z) 



3 
V2 



(13) 



which may be derived from the full Hamiltonian Eq. ([8]l and captures many essential 
elements of the model; for example, the dispersion Eq. ( fTOl i can be found immediately 
from det (-ffeff(^) ~ z) = Q (see Ref. [25|). The process of deriving this effective 
Hamtilonian Eq. ( fTSI ) from the full Hamiltonian for our prototype model is outlined 
in App. [Aland represented diagrammatically in the transition from Fig.lTfa) to (b). 

From here, one might naively expect that the number of EPs in OQS is now deter- 
mined by the dimension rie^s of H^a according to r := rieff (rieff — 1) and indeed this 
(by chance) gives the correct answer for our prototype model from rieff — 2. How- 
ever, this is misleading; one must emphasize that iJeff itself contains z-dependence, 
which results in a characteristic equation of the form H^fi{z)ipz = z^^ that is usually 
of higher order than the dimension of H^ff itself. In fact, it can be shown that in the 
case of an effective Hamiltonian i/etr of dimension n with quadratic dispersion in the 
leads, the number of eigenvalues is given by 2?! and that the present model is just a 
special case in which two of the solutions are trivial (see App. [Bjof the present text 
as well as Apps. D and H of 11241 ). 

Indeed, observing the form of Eq. (fT3] l it is straight-forward to realize that the 
dimension of iJgff increases by 1 for each discrete impurity we add to the system or 
for each distinct chain (continuum) we add. However, we must take into account that 
introducing a new chain also doubles the number of polynomial eigenvalue equations 
that must be solved, effectively doubling the number of solutions. Hence the number 
of solutions in the general case is given by 2""^ (no + nc) in which rio is the number 
of discrete impurity sites in Hy) and nc is the number of (non-degenerate) continuum 
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channel^ in Hq- Therefore, the number of EPs in OQS with quadratic dispersion is 
given by 

r = 2"^(nD + nc)[2"^K + nc) - 1]. (14) 

For easy-to-check examples consider the systems in Refs. II17II21||221 . 

Returning now to our present evaluation, Eq. ( fTOl i can be immediately solved for 
the two eigenvalues 



' 1 - 2g2 ^ 1-252 



z±[ed,g) = ed- —^±g- —^ . (15) 



We write the associated wave functions along the chain sites in the (unnormalized) 
form tp±{x) = {ip±\ x) ~ e*'^±^ 123)1 with the effective /c-value defined as k± = 
cos~^(— z±) and given explicitly as 



i^^e,±^4-il-2g^)^ 



A:±=ilog(-ed±Vq-(l-V)l- (16) 

Note that at the EPs we have fc+ — fc_ such that the two eigenfunctions collapse into 
one another. 

Focusing on the EPs ej = e±, we briefly illustrate the method outlined for an- 
alyzing EPs in OQS by plugging S{z) from Eq. ( fTOl i into the condition Eq. (|6]l to 

write 

as a condition to find the two values for the center, with q ^ ± a sign designation to 
be fixed below. Equation (fTTT i easily yields 

z^ = ± , -^ ■ (18) 

We perform a consistency check to find that the sign q = + (q = —) before the root 
in Eq. (fTTT i is associated with z~ (z+). Making consistent use of this sign conven- 
tion and plugging the value z = z^ into the original dispersion Eq. ( fTOl i yields the 
condition for the EPs in parameter space as ej = e±, just as we found in Eq. (fT2l i 
from /2(ed, g) = 0. Finally we use Eq. (fTSl l to write the center values in terms of the 
other system parameter ej as z^ ~ ±(1 + ed2)/2 ej. From here we could continue 
the method outlined in the text following Eq. (|6]l to obtain the eigenvalue expansion 
Eq. (HJi. For a model with a more complicated self-energy function X!{z) this would 
probably be necessary but in the present case we may obtain the full expansion more 
easily. 



4 Prototype spectrum and eigenvalue expansion 

Before obtaining the fractional power series, we discuss the spectrum in the case 
g < l/\/2 and demonstrate that the resonant state emerges precisely at the EPs 



^ For an example of a multi-channel model, consider the two-chain model in 1171 . in which it is straight- 
forward to derive the effective Hamiltonian H^ff of dimension 3, giving rise to 12 solutions. 
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Fig. 2 (Color online) (a) Real part of the 2-j- spectrum (both eigenvalues), (b) imaginary part of 2_ 
(proportional to decay rate), (c) real part of the effective wave number A;_ for the 2_ eigenvalue, and (d) 
imaginary part of the effective wave numbers k± for both eigenvalues, plotted in each case as a function 
of Ed for fixed g = 0.67 < 1/V2. 



e±. The real and imaginary parts of the spectrum are plotted in Figs. |2a) and (b), 
respectively, as a function of ej for fixed g = .67. Solutions in the first (second) 
Riemann sheet are shown with a solid (dashed) curve. In Fig. |2ta) we find for the 
case Ed > e^ with 



±(1 



(19) 



there exists exactly one bound state given by the purely real solution z_ > 1; this 
bound state occurs in the first Riemann sheet in z space. For this case we find there 
also exists a real solution in the second Riemann sheet that we call an anti-bound 
state; we discuss the properties of such a state momentarily. If we decrease ea to 
the value ej = e^ (> e+ > 0) the bound state in the first sheet turns around the 
branch point z = 1 (i.e., it brushes against the energy continuum) before crossing 
into the second Riemann sheet as well Ill6ll29|[30l . The crossover point is indicated 
in Fig. |2ja) with a dotted brown circle. For the negative side of the spectrum a similar 
process occurs in the vicinity of e^ — e^ {< e^ < 0), in which z_ and z+ switch 
roles. 

This transition can be understood in terms of the evolution of the imaginary part 
of fc_, which from Eq. (fTSI l is easily shown to be positive-valued for ea > e^, as 
shown in Fig.|2d). This means that the wave function ip- (x) = e''^^^ is bounded for 
all X (the electron tends to occupy levels close to the impurity). If we decrease ej the 
imaginary part of fc_ vanishes at ed = ej^ and then becomes negative for e^ < e^, and 
hence the wave function e**^-^ has become unbounded. Clearly, this type of solution 
no longer yields an ordinary bound state and hence we call it an anti-bound state (or 
virtual state lfT6l ) in the second Riemann sheet in z space; in this 'state' the electron 
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tends to occupy levels along the chain far away from the impurity. The properties of 
the present model that induce this transition are detailed in Ref. 1^). 

Physically, the transition from a bound to an anti -bound state is associated with an 
amplification of the long-time deviations from exponential decay 1291 that are well 
known to occur in quantum systems 01ll32|[33l . This transition is also associated 
with an amplification in the local density of states profile at the impurity site Il34ll . 

We emphasize that e^ do not represent EPs since there is no coalescence of the 
wave functions here; in fact, these points are more closely related to the spectral sin- 
gularities II15II36I that may occur in the continuous spectrum of OQS if we relax our 
conditions on the Hamiltonian to allow for explicit non-Hermiticity. From a physical 
perspective, it has been shown that a diverging or vanishing reflection coefficient is 
obtained precisely at the point in parameter space at which the spectral singularity 
occurs. These special points may be found in our prototype model by restricting the 
impurity energy e^ to complex values. For a detailed analysis of the spectral singular- 
ities occurring in the present model, consult Ref. (T5\ . 

Now we discuss the emergence of the resonance state at the real-valued EPs e^ = 
e± (highlighted with solid green circles in Fig. |2]i, which appears to be a general 
phenomena. As shown above, in the region of the spectrum e+ < ed < el^ there are 
two anti-bound states in the second Riemann sheet in z space. At the precise values 
Ed = e+ a different transition occurs as these two eigenvalues coalesce giving rise to 
a complex conjugate pair of solutions; this fact is apparent from the root appearing 
in the eigenvalue functions Eq. ( fTSl l. The eigenvalue z_ gives the resonant state with 
negative imaginary part (decay width), plotted in Fig. lUb). Hence we find that the 
emergence of the resonance (associated with time irreversibility) occurs precisely at 
the EP We finally point out in Fig.|2jc) that for all ed < e^ (ej > e^) the real part of 
k^ is (tt); for the anti-bound states, the wave vector always lies at one edge of the 
Brillouin zone. It is only when the resonant state appears that fc_ takes a value inside 
the unperturbed channel of fc-states, k G (0, tt); from Eq. (fT6b we find the exact value 
is given by k^es = ikc + 4>ies with kc = log ■\/l — 2g2 and 

arctan(|A^/^|/ed) , e_ < ed < 0, 

(20) 
arctaii(-|Ai/2|/ed) +7r, < ed < e+, 

in which A = ed^ —(1 — 2g^) (see Figs. |2jc,d)). Finally, for completeness we plot 
Im 2:+ and Re fc_|- in Figs.[3a,b), respectively. 

We would now like to demonstrate the process of adiabatically encircling the EP. 
Up until this point we have considered the chemical potential ed at the impurity site 
to take on strictly real values. However, let us for a brief moment consider the pos- 
sibility of a complex- valued Cd, as has been studied for the present model in Ref. 
ifTSl . (We emphasize that in this reference the goal of the author is, unlike the present 
situation, to consider a complex chemical potential from the start, and then investi- 
gate the consequences of that generalization.) We further imagine that the chemical 
potential at the impurity may be controlled by a gate potential, and hence we vary 
the impurity energy Cd around the EPs ± -^/l — 2g'^ for a fixed g value according to 
(ed(0))^ = 1 - 2g'^ + (5e*^, with S < 2{1 - 2g^) to avoid encircling the second 
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Fig. 3 (Color online) (a) Imaginai'y part of 2+ and (b) real part of the effective wave number fc+ , plotted 
as a function of t^ for fixed g = 0.67 < l/\/2. 



EP as well. Applying this parameterization to the eigenvalues Eq. ( fTSl ) we see that 

V (ed(^)) — (1 — 29^) — \/je*^/^, which illustrates that z_ — > z^ and vice versa 
as 9 is varied from to 2tt. Hence the cycle structure Eq. (|2]i is given by {z_ ,2:+} in 
the vicinity of either EP for the present model. 

In order to demonstrate the utility of our formalism and to motivate the succeed- 
ing QPT analogy we now obtain the fractional power eigenvalue expansion Eq. (|4| 
following a method from Ref. J37J. To do this we encode the variation of the system 
parameters about the EPs, in this case using ^^ = A = e^^ —(1 — 2g^); we aim to 
perform expansions about the non-analytic points ej = e±, but the non-analyticity 
prevents us from doing so directly. To overcome this problem, we write a new func- 
tion for both eigenvalues s = ± at either EP as 



0s(e±,C) = ^s(ed = e±,ff), 



(21) 



in which we allow g ~ J Q"^ + 1 — ej_/\/2 to vary as a function of C,\ this effec- 
tively measures the distance from the EP in parameter space. We can expand this 
new function 0s in an ordinary Taylor series in C, and then use the reverse parameter 
transformation to write the original eigenenergies as 



Zs{ei,g) = z^ + s 



1 ?2 
i-^Ai/2 
24 



1 

2il 



00 

E 

n=2 






(22) 



with A = /2(ed, .9) for comparison with Eq. dUi. Keep in mind that each occurrence 
of l± above specifies one of the two EPs, while s = ± refers to one of the two 
eigenvalues (i.e., Eq. ( l22l i is actually four equations). This expression provides the 
realization of the generic expansion Eq. (|4|i for the case of the semi-infinite chain 
with an endpoint impurity in which case p takes the value p = 2. 



5 QPT analogy for the appearance of the resonance at the EP 



As a final comment on the EPs in OQS, we hypothesize a QPT analogy for the emer- 
gence of irreversibility at the EPs where the resonance appears. This analogy is sim- 
ilar to that proposed in Ref. lfT9l . in which the decay rate F = — 2 Im z_ is related to 



12 Savannah Gai'mon et al. 

the order parameter. We use the EP e+ as our example; the analogy works out sim- 
ilarly for the negative side of the spectrum involving the EP e_ . Let us focus on the 
'non-critical case' e+ < ed < e^, where z± are both real. Here we have a pseudo- 
energy gap between the two anti-bound states given by Az = z+ — z^ ^ A^/^ ^ 
(ed — e+)^/^ as shown in Eqs. (fTST i or ( l22l i: this pseudo-gap vanishes at the 'critical 
point' Cd = e+, see Fig. |2a). Then below the transition e^ < e+ the decay width 
appears as rj-es ~ A^/^ ~ (ed — e+)^/^, representing broken time-reversal symmetry. 
We may carry the analogy further by defining a function that describes the correla- 
tions between the impurity site and an arbitrary chain site x through the resonant state 
by 

CUx-d) = ^.(f) dz{x\—^\d), (23) 

where C_ represents a contour in the complex z plane that surrounds only the eigen- 
value z_ lim . We evaluate this to find 



C,-es(x; d) - e*''-^ = e'°s Vi-V^e"^"»^. (24) 

We recall that on the 'non-critical' side ed > e+, -Z- is an anti-bound state with 
4> = IT, giving purely real-valued correlations. Meanwhile, on the 'critical' side < 
Ed < e+ we have (fires ~ arctan (— |A^/^|/ ed) + ti" from Eq. ( l20b . (plotted in Fig. 
|2|c)), which gives complex-valued correlations. Further, we see that 0ies takes a value 
within the continuum of fc-states 0,-es G (0, tt) (the diagonalized states for the system 
if no impurity were present), signifying that the resonance state escapes the impurity 
and propagates along the chain with decay width ries. 

This is reminiscent of the dynamical phase transition in Ref. 1191 that occurs 
when a resonance state, coupled to a large number of otherwise similar resonance 
states, suddenly transitions from a narrow resonance to a wide resonance. The order 
parameter associated with this process is proportional to the resonance decay width, 
which suffers a discontinuity in the first derivative at the transition point. Simultane- 
ous with this transition, the wave function of this state "collects all the corresponding 
components from the wave functions of all the other states." 

In the present case we instead have two anti-bound states (zero decay width) 
that spontaneously give rise at the EP to a resonance (finite decay width) and anti- 
resonance pair. Meanwhile the correlations for the resonance appear similar to what 
we would find for the unperturbed k states in the uncoupled model 

Co{x; k) = {x\ 77— IT 1^^) ~ ^''"' (25) 

z — Uu — tic 

with k G (0, tt). Hence the resonance state in some sense mimics the |fc) states, which 
are themselves nothing more than a collective oscillation of all the \x) states in the 
original site basis. From Eq. (l24l i we now introduce an inverse correlation length on 
the critical side of the transition as ^^^ = 0res; further we note ^ ~ tt + X~^'^ ~ 
TT + (cd — e-|_)^^/'^. Since the decay width in this case also suffers a discontinuity 
in the first derivate at the EP, we interpret this as the order parameter varying as 
r ~ ^~^ ~ (ed —e^y^^, which yields the dynamic exponent z = 1 1381 : we propose 
that this exponent might be interpreted as representing the time dimension in which 
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the symmetry is spontaneously broken with the appearance of exponential decay in 
association with the resonance state. 

We note our argument above does not rest on the details of the prototype model 
but rather on generic features associated with the complex eigenvalue; for example 
a similar transition occurs for a quite different physical situation in Ref. IS] and our 
QPT analogy should directly apply in the systems studied in Refs. II14II17||T1II22||241 
[34l . Indeed it seems that the resonant state in OQS should always appear at a real- 
valued EP described by the expansion in Eq. (|4|i. 

Recently, an experimental demonstration of precisely this type of transition has 
been performed in an electronic system consisting of a coupled LRC circuit pair ll39ll . 
In that system, one LRC circuit is designed for amplification while the other provides 
a compensatory gain, resulting in a system with manifest parity-time T'T-symmetry 
Pol . By varying the system parameters, the authors of Ref. ||39l are able to demon- 
strate a transition from real-valued eigenfrequencies (the "exact T'T-symmetric phase") 
to complex eigenfrequencies (the "broken T'T-symmetric phase"). 

It is worth noting at this point that although the Hamiltonian for our prototype 
model given in Eq. (O appears to be Hermitian, there is however an underlying non- 
Hermiticity associated with the continuum component of the system. On a deeper 
level the system is in fact T'T-symmetric rather than Hermitian. This underlying 
T'T-symmetric nature is revealed for the present system by writing the effective 
Hamiltonian iJeff in Eq. ( fT3] l (also see iJgff re-written in terms of the k variable in 
Eq. (ISTTi). Indeed, applying the Hermitian transformation on Eq. dSTT i (with F = 1/2 
for the prototype model) reveals that the system is non-Hermitian [251, and indeed the 
non-Hermiticity comes from the lower-right entry for iJgff, which is the self-energy 
associated with the continuum. However, ifeff is indeed symmetric under the VT 
transformation. This underlying T'T-symmetric property is again a generic feature of 
OQS. 



6 Conclusion 

In this paper we have presented the basic mathematical properties of exceptional 
points occurring in the spectrum of open quantum systems, along with a technique for 
their location and analysis in these systems. This technique includes our method for 
locating EPs in OQS detailed in Sec.|2] While the location method we have presented 
is oriented towards analytically solvable models, it nevertheless may prove useful in 
some numerical situations assuming that we may at least explicitly write the self- 
energy function S{z) in integral form. We have also outlined our method to obtain 
the fractional power expansion given in generic form in Eq. (JUi, which itself has been 
generalized from the expression appearing in Kato's text Ref. f\3\ (apparently this 
generalization is required to accommodate OQS, if perhaps not wider circumstances). 
We have further applied these techniques to the case of a semi-infinite chain with 
an endpoint impurity. We found that the resonant state eigenvalue emerges precisely 
at a real-valued EP, which should be a general result for OQS. Further we offered 
a quantum phase transition analogy for the process in which one crosses this real- 
valued EP in parameter space. According to this analogy the decay width, which ex- 
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periences a discontinuity in its first derivative at the EP, plays the role of order param- 
eter As the resonant state gives rise to time irreversibility via exponential decay, this 
appears to be related to the time-symmetry breaking mechanism in PT-symmetric 
systems, including OQS. 

A Derivation of the effective Hamiltonian Hett for the prototype model 

In order to keep our presentation reasonably self-contained, we briefly outline the method from Ref. [251 
to obtain the effective Hamiltonian reported in Eq. U3) from the full Hamiltonian for our prototype model 
given in Eq. (8). We start by writing the wavefunction along the chain in generic form as 

tp{x) = Ce^''^ for a: > 2 (26) 

with C an undetermined normalization constant. To solve the model from this point we find it useful to 
close the Schrodinger equation on the right with an arbitrary bra {x\ to obtain 

(x\H\iP) = e(x\^) (27) 

in which H is the full Hamiltonian from Eq. (8), while x might represent the impurity site or any site 
along the chain x £ {d, 1, 2, . . .}. If we specify any site x > 2 and plug Eq. {8) into Eq. )27t then, after 
applying the usual commutation relations, we find 

- i ( (x - l\i,) + {x + l|t/.> ) = ekixm. (28) 

Then plugging in Eq. J26) for the wave function we immediately find e^ = — cos k for the dispersion 
along the chain. Next we solve for the normalization constant C by choosing a; = 2 in Eq. J27K which 
yields 

-i((l|V'>+Ce3''^)=efeC7e2*. (29) 

After applying efc = — cos k here we re-arrange to find C = {l\ip)e~'''. 

Finally, evaluating Eq. (27) for the cases x = 1 and x = d we obtain the matrix equation 



« ti)(5?i!;n=4 5f!!n, (30) 






in which e = z is the discrete eigenvalue associated with the impurity sector. Applying the transformation 
k = cos~^(— z) in the 2x2 matrix on the LHS of this equation we find — e'*' = 2 — Vz"^ — 1 = -^(2), 
hence we obtain the effective Hamiltonian H^ff reported in Eq. )13) in the main text. We now obtain the 
characteristic equation by solving dot (z — Hci[{z)) = 0, as reported in Eq. UOt . 



B Comment on the number of discrete solutions resulting from a general 
effective Hamiltonian Hett 

In this Appendix we briefly comment on the number of solutions emerging from the characteristic equation 
under the effective Hamiltonian treatment (the reader should also consult App. D and H of 1241 for similar 
discussions). Let us consider a slightly more general form for H^-ff of dimension two as 

HM^) = ( _\ _^,k ]■ (31) 

(compare this with Hcff from Eq. BO) for the semi-infinite chain model with an endpoint impurity). 
Note that we may just as well solve the characteristic equation in terms of k rather than z by writing 
dot (z{k) — Hcff(k)) = 0, which gives 

(e''= -I- e-''' -I- 2 ed) ((1 - 2F)e* -|- e~*) - 2g'^ = (32) 
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after applying the transformation z(fe) = — cosfc = (e*'°+e~''°)/2 foratight-binding chain. Making the 
replacement to = e**^ here and multiplying through by w^ immediately yields a fourth order polynomial 
equation for to, hence a general effective Hamiltonian of dimension two should be expected to yield four 
discrete eigenvalues. From this point it is not difficult to convince oneself that adding a single discrete 
impurity site increases the dimension of H^ff by one and simultaneously introduces two new solutions to 
the characteristic equation. Meanwhile, adding a new channel also increases the dimension of H^ff by one 
but also requires us to solve a second polynomial equation, hence doubling the total number of solutions 
after the increase in the dimension of the effective Hamiltonian has already been taken into account. 

However, our semi-infinite chain model with an endpoint impurity happens to have the special value 
F = 1/2, which we note is the only value for F under which the order of the polynomial Eq. (32) is 
reduced from four to two. Hence our present model represents a special case with regard to the number of 
discrete solutions (and therefore the number of EPs as well). 



Acknowledgements We would like to thank E. C. G. Sudarshan, G. Ordonez, J. Pixley and S. Kirchner 
for stimulating discussion on the QPT analogy and T. Petrosky for helpful comments on an earlier draft. 
The research of S. G. was supported by CQIQC, the Sloan research fellowship of D. S. and MPI-PKS. The 
research of N. H. was supported by Grant-in-Aid for Scientific Research No. 173401 15 from the Ministry 
of Education, Culture, Sports, Science and Technology of Japan. 



References 

1. W. D. Heiss: Phases of wave functions and level repulsion. Eur. Phys. J. D 7, 1 (1999). 

2. C. Dembowski, H.-D. Graf, H. L. Harney, A. Heine, W. D. Heiss, H. Rehfeld and A. Richter: Experi- 
mental observation of the topological structure of exceptional points. Phys. Rev. Lett. 86, 787 (2001) 

3. C. Dembowski, B. Dietz, H.-D. Graf, H. L. Harney, A. Heine, W. D. Heiss, and A. Richter: Observation 
of a chiral state in a microwave cavity. Phys. Rev. Lett. 90, 034101 (2003). 

4. S.-B. Lee, J. Yang, S. Moon, S.-Y. Lee, J.-B. Shim, S. W. Kim, J.-H. Lee, and K. An: Observation of 
an exceptional point in a chaotic optical microcavity. Phys. Rev. Lett. 103, 134101 (2009). 

5. I. Rotter: A non-Hermitian Hamilton operator and the physics of open quantum systems. J. Phys. A: 
Math. Theor. 42. 153001 (2009). 

6. J. Rubinstein, P. Sternberg, and Q. Ma: Bifurcation Diagram and Pattern Formation of Phase Slip 
Centers in Superconducting Wires Driven with Electric Currents, Phys. Rev. Lett. 99, 167003 (2007). 

7. N. Moiseyev, and S. Friedland: Association of resonance states with the incomplete spectrum of finite 
complex-scaled Hamiltonian matrices, Phys. Rev. A 22, 618 (1980). 

8. R. Lefebvre, O Atabek, M. Sindelka, and N. Moiseyev: Resonance Coalescence in Molecular Photodis- 
sociation, Phys. Rev. Lett. 103, 123003 (2009). 

9. H. Cartarius, J. Main, and G. Wunner: Exceptional points in atomic spectra, Phys. Rev. Lett 99, 173003 
(2007). 

10. E. Hernandez, A. Jauegui and A. Mondragon: Non-Hermitian degeneracy of two unbound states, J. 
Phys. A: Math. Gen. 39, 10087 (2006). 

11. M. R. Zimbauer, J. J. M. Verbaarschot and H. A. Weidenmilller: Destruction of order in nuclear 
spectra by a residual GOE interaction, Nucl. Phys. A 411, 161 (1983). 

12. P. Cejnar, S. Heinze and M. Macek: Coulomb Analogy for Non-Hermitian Degeneracies near Quan- 
tum Phase Transitions, Phys. Rev. Lett. 99, 100601 (2007). 

13. T. Kato, Perturbation Theory for Linear Operators, Springer- Verlag, Berlin (1980), pp. 62-66. 

14. T. Petrosky, C.-O. Ting, and S. Garmon: Strongly Coupled Matter Field and Nonanalytic Decay Rate 
of Dipole Molecules in a Waveguide, Phys. Rev. Lett. 94, 043601 (2005). 

15. S. Longhi: Spectral singularities in a non-Hermitian Friedrichs-Fano-Anderson model, Phys. Rev. B 
80, 165125 (2009). 

16. R. A. Bustos-Marun. E. A. Coronado, and H. M. Pastawski: Buffering plasmons in nanoparticle 
waveguides at the virtual-localized transition, Phys. Rev. B, 82, 035434 (2010). 

17. S. Garmon, H. Nakamura, N. Hatano, and T. Petrosky: Two-channel quantum wire with an adatom im- 
purity: Role of the van Hove singularity in the quasibound state in continuum, decay rate amphfication, 
and the Fano effect Phys. Rev. B 80, 1 15318 (2009). 

18. E. N. Bulgakov, 1. Rotter, and A. F. Sadreev: Phase rigidity and avoided level crossings in the complex 
energy plane Phys. Rev. E 74, 056204 (2006). 



16 Savannah Gai'mon et al. 



19. C. Jung, M. Muller, and I. Rotter: Phase transitions in open quantum systems, Phys. Rev. E 60, 114 
(1999). 

20. S. Longhi: Nonexponential decay via tunneling in tight-binding lattices and the optical Zeno effect, 
Phys. Rev. Lett. 97, 1 10402 (2006). 

21. S. Longhi: Boimd states in the continuum in a single-level Fano- Anderson model, Eur Phys. J. B 57, 
45 (2007). 

22. S. Tanaka, S. Garmon, G. Ordonez, and T. Petrosky: Electron trapping in a one-dimensional semicon- 
ductor quantum wire with multiple impurities, Phys. Rev. B 76, 153308 (2007). 

23. N. Hatano, K. Sasada, H. Nakamura and T. Petrosky: Some Properties of the Resonant State in Quan- 
tum Mechanics and Its Computation, Prog. Theor Phys. 119, 187 (2008). 

24. K. Sasada, N. Hatano, and G. Ordonez: Resonant spectrum analysis of the conductance of an open 
quantum system and three types of Fano parameter, J. Phys. Soc. Japan 80, 104707 (2011). 

25. N. Hatano: Equivalence of the effective Hamiltonian approach and the Siegert boundary condition for 
resonant states, to be pubHshed in Fortschritte der Physik, DOl: 10.1002/prop.201200064 (2012). 

26. E. Jarlebring, S. Kvaal, and W. Michiels: Computing all pairs (A, /i) such that A is a double eigenvalue 
of A + IJ.B KU Leuven. TW559 (2010). 

27. R. Lefebvre andN. Moiseyev: LocaUzation of exceptional points with Pade approximants, J. Phys. B: 
At. Mol. Opt. Phys. 43, 095401 (2010). 

28. R. Uzdin and R. Lefebvre: Finding and pinpointing exceptional points of an open quantum system, J. 
Phys. B: At. Mol. Opt. Phys. 43, 235004 (2010). 

29. S. Garmon, T. Petrosky, Y. Nikulina, and D. Segal: Amplification of non-Markovian decay due to 
bound state absorption into continuum, to be published in Fortschritte der Physik, arXiv:1204.6141 

30. E. N. Economou: Green's Functions in Quantum Physics, 2nd ed., Springer- Verlag, Berlin (1983). 

31. C. B. Chiu, B. Misra, and E. C. G. Sudarshan: The evolution of unstable states and a resolution of 
Zeno's paradox, Phys. Rev. D 16, 520 (1977). 

32. L. Fonda, G. C. Ghirardi, and A. Rimini: Decay theory of unstable quantum systems. Rep. Prog. Phys. 
41, 587 (1978). 

33. J. Martorell, J. G. Muga, and D. W. L. Spring: Quantum post-exponential decay, Lect. Notes Phys. 
789, 239 (2009). 

34. A. D. Dente, R. A. Bustos-Mariin, and H. M. Pastawski: Dynamical regimes of a quantum SWAP gate 
beyond the Fermi golden rule. Phys. Rev. A 78, 0621 16 (2008). 

35. B. A. Fuchs and V. I. Levin: Functions of a Complex Variable, Vol. 11, Pergamon Press (1961). 

36. A. Mostafazadeh: Spectral Singularities of Complex Scattering Potentials and Infinite Reflection and 
Transmission Coefficients at Real Energies, Phys. Rev. Lett 102, 220402 (2009). 

37. K. Knopp: Theory of Functions Pt. 2, Dover Publications, New York, pp. 126-129. 

38. S. Sachdev, Quantum Phase Transitions, Cambridge University Press (2001). 

39. Joseph Schindler, Ang Li, Mei C. Zheng, F. M. Ellis, and Tsampikos Kottos: Experimental study of 
active LRC circuits with "PT-symmetries, Phys. Rev. A 84, 040101(R) (201 1). 

40. C. M. Bender and S. Boettcher: Real spectra in non-Hermitian hamiltonians having T'T-symmetry 
Phys. Rev. Lett. 80, 5243 (1998). 



